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ABSTRACT 

ISO-Oph-50 is a young low-mass object in the ~ 1 Myr old Ophiuchus star forming 
region undergoing dramatic changes in its optical/near/mid-infrared brightness by 2- 
4 mag. We present new multi-band photometry and near-infrared spectra, combined 
with a synopsis of the existing literature data. Based on the spectroscopy, the source 
is confirmed as a mid M dwarf, with evidence for ongoing accretion. The near-infrared 
lightcurves show large-scale variations, with 2-4 mag amplitude in the bands IJHK, 
with the object generally being bluer when faint. Near its brightest state, the ob¬ 
ject shows colour changes consistent with variable extinction of A Ay ~ 7 mag. High- 
cadence monitoring at 3.6 pm reveals quasi-periodic variations with a typical timescale 
of 1-2 weeks. The best explanation for these characteristics is a low-mass star seen 
through circumstellar matter, whose complex variability is caused by changing inho¬ 
mogeneities in the inner parts of the disk. When faint, the direct stellar emission is 
blocked; the near-infrared radiation is dominated by scattered light. When bright, the 
emission is consistent with a photosphere strongly reddened by circumstellar dust. 
Based on the available constraints, the inhomogeneities have to be located at or be¬ 
yond ~ 0.1 AU distance from the star. If this scenario turns out to be correct, a major 
portion of the inner disk has to be clumpy, structured, and/or in turmoil. In its ob¬ 
servational characteristics, this object resembles other types of YSOs with variability 
caused in the inner disk. Compared to other objects, however, ISO-Oph-50 is clearly 
an extreme case, given the large amplitude of the brightness and colour changes com¬ 
bined with the erratic behaviour. ISO-Oph-50 has been near its brightest state since 
2013; further monitoring is highly encouraged. 

Key words: stars: low-mass, brown dwarfs; stars: activity; stars: pre-main-sequence; 
accretion, accretion discs 


1 INTRODUCTION 

Variability has long been known as a characteristic ob¬ 
servational signature of young stars. Most common are 
lightcurves with periods of days to weeks due to spots on 
the stellar surface, a fact that has been used for almost 
three decades to infer the rotation periods of young stars 


(Bouvier & Bertout 1989 Herbst et al. 20071. Apart from 


these regular changes, a fraction of young stellar objects 
(YSOs) shows evidence for large-scale, partly irregular vari¬ 
ations (Herbst et al. 19941. This was originally one of the 


criteria for the definition of the class of ’T Tauri stars’ (Joy 


19451. Some of the irregular variability follows patterns that 


E-mail: asllO@st-andrews.ac.uk 


are well-explained by astrophysical phenomena, including 


hot accretion-induced spots (Vrba et al. 1989 

Bouvier et al. 

19951, magnetic flares (Rydgren & Vrba 1983 

Welty & Ram- 


sey|1995 1, accretion bursts (the FU Ori and EX Lupi events; 


Herbig||1989; Hartmann & Kenyon 1996), and obscurations 


by circumstellar material. The latter includes the UX Ori 
phenomena for Herbig Ae/Be stars ( Grinin et al. 1 199T ), but 
also other types (Hamilton et al.||2001 Rodriguez-Ledesma 


et al. 2012). Variable extinction may be the reason for the ex¬ 


treme fluctuations observed in the prototype, T Tauri (Beck 
& Simon 2001). Often young stars are affected by several of 


these processes. For a comprehensive discussion of the vari¬ 
ous factors that contribute to YSO variability and a discus- 


et al. (2009); Wolk et al. (2013). 


sion of possible models, see Carpenter et al. (2001); Scholz 
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Although extreme variables have been defining proto¬ 
types for young stellar objects, their overall fraction is low. 
In the near-infrared, only 2-3% of an unbiased sample of 
YSOs show variations larger than 0.5 mag over timescales 
up to several years ( Scholz||2012 1. In the optical, that frac¬ 
tion rises to ~ 10% (Rigon et al., in prep.). The subset of 
YSOs variable on timescales of years is only 10-20% in the 
mid-infrared (Rebull et aLj20141. The fraction of extremely 
variable YSOs drops with age; only very few arc known older 
than 5Myr ( Mamajek et al.|2012 Grinin et al.|2010| , in line 
with typical disk lifetimes. Despite their rarity, YSOs with 
strong variations present unique insights into specific pro¬ 
cesses and phases occuring during the first few million years 
in a star’s life (e.g., episodes of strong accretion). In partic¬ 
ular, YSOs which are obscured by their disks can provide 
information about the structure and the evolution of cir- 
cumstellar accretion disks, the properties of the dust grains 
in the vicinity of the star, as well as the initial conditions of 
planet formation (e.g. Mora et al.||2002' Herbst et al. ||2008 


Rodriguez et al. |2013|). 

So far, most of the detailed studies of extreme T Tauri 
variability have been carried out in the optical and for stars 
more massive or similar in mass to the Sun. Current stud¬ 
ies are driven by large-scale optical variability surveys (e.g., 
Super-WASP, KELT, PTF, Kepler, Gaia). Many deeply em¬ 
bedded and very low mass YSOs are not covered by such 
surveys. On the other hand, given the fact that most stars 
have masses < 0.5 Mg and that the occurence of extreme 
variability declines with age, it is interesting to investigate 
the aforementioned phenomena for lower-mass objects and 
in the embedded phase, to probe how applicable the pro¬ 
posed explanations are to the youngest and more typical 
YSOs. As an example, it was recently found that an embed¬ 
ded very low mass star undergoes an FU Ori-likc outburst, 
demonstrating that accretion bursts can occur over a wide 
range of stellar masses ( Caratti o Garatti et al.|201 1). 

Here we present new observations for the infrared source 
ISO-Oph-50, located in the ~ lMyr p Op h star form i ng re- 
gior0 The object was identified by Bontemps et al. (20011 
from ISO observations and tentatively classified as Class III 
(i.e. without disk) based on the infrared SED slope. More 
recent studies based on Spitzer data put it at an earlier evo¬ 
lutionary state of Class I or II (|Gunther et al.|[2014| |Alves 


de Oliveira et al. 2012). ISO-Oph-50 has been known to 


show extreme variability in the infrared (Alves de Oliveira 
|fe Casali|2008| [Scholz et al.|2()l 3~j ). So far, the only available 
spectrum shows that the object is an M-dwarf (Alves de 
Oliveira et al.|[2012). The source has excess emission in the 


mid-infrared indicating the presence of circumstellar dust, 
but so far there is no evidence for ongoing accretion. In 
this paper we analyse new multi-band photometry and near- 
infrared spectra and discuss possible explanations for the 
variability. 


1 2MASS J16263682-2419002, J2000 coordinates: 16:26:36.8, - 
24:19:00 


2 OBSERVATIONS 


2.1 SMARTS photometry 


We used the 1.3 m telescope at the Cerro Tololo Interna¬ 
tional Observatory to obtain long-term photometry for ISO- 
Oph-50, in the framework of guaranteed observing time 
through the SMARTS collaboration. In total we observed 
the source in 20 nights spread over 2 years. The telescope 
is equipped with Andicam, a dual-channel camera which 
takes optical and near-infrared images simultaneously. Our 
observing scripts first take images in I- and J-band, then 
R- and K-band (with the exception of the very first epoch 
in which we observed in H instead of K). ISO-Oph-50 is 
not visible in the R-band frames, and only barely detected 
in most I-band images. All observations were done at low 
airmass ranging from 1.0 to 1.5. 

In the near-infrared bands, the object was observed in a 
5 position dither pattern with integration times of 5 x 30 sec 
per band. A standard image reduction was carried out, in¬ 
cluding sky subtraction, flatfielding and aperture photom¬ 
etry. While the J-band images show only our target, the 
H- and K-band images cover three more stars with 2MASS 
photometry. We calibrated the photometry in relation to the 
2MASS values for these three stars. The offset between in¬ 
strumental K-band magnitudes ivi na t and 2MASS K-band 
magnitudes Kim remains stable at Ai ns t — A 2 M = 3.0 ± 0.1 
and does not significantly depend on the colour of the com¬ 
parison stars. This offset includes the constant zeropoint 
and an extinction term, the variations of ±0.1 are consis¬ 
tent with a combination of photometric errors and the vary¬ 
ing airmass. Typical errors of the K-band magnitudes are 
±0.05 mag, combining the calibration uncertainties and the 
photometric error. 


Since we do not have a calibration star in the J-band, 
we have established the offset between instrumental mag¬ 
nitudes and 2MASS magnitudes from observations of other 
fields, published in Scholz et al. (2013). The offset in the 
J-band is Ji ns t — J 2 M ~ 2.5 ±0.1 through seasons 2012 to 
2014. Similarly to the K-band, this value is stable within the 
quoted errorbars and for the given airmass range. We adopt 
this offset for all our observations of ISO-Oph-50. Since we 
cannot measure the offset for the images themselves, the re¬ 
sulting uncertainties may be as large as ±0.1 (comparable to 
the variations in the offset). The calibrated J- and K-band 
lightcurves are plotted in Fig. [T] 


For the optical images, we obtained three images per 
band at the same position, each with an integration time of 
110 sec. To increase signal-to-noise, the I-band images were 
co-added. As calibration star, we used the object 2MASS 
J16264285-2420299 (or ISO-Oph-62), which has a DENIS 
I-band magnitude of 14.355 and does not show signs of sig¬ 
nificant variability. The uncertainty in the I-band measure¬ 
ments is in the range of ±0.3 mag. This estimate includes 
the photometric error (±0.1 mag) plus possible calibration 
systematics. 


The full SMARTS photometry is listed in the appendix 
in table |A1[ This table also includes the first epoch of 
SMARTS photometry in J- and H-band which was already 
published in Scholz et al. (20131 for the 12th of August 2012. 
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Figure 1. SMARTS J/K-band lightcurve for ISO-Oph-50 cov¬ 
ering two years. Errors in J and K are ±0.1 and ±0.05 mag, re¬ 
spectively. 



Figure 2. NTT/SOFI spectra for ISO-Oph-50 taken on April 
13th and 19th 2014. The two spectra are calibrated against each 
other using the K-band magnitudes obtained from the acquisition 
images. The accretion-related lines Paschen 0 and Bracket 7 are 
marked. 


2.2 NTT/SOFI Spectroscopy 

We took two spectra of ISO-Oph-50 using SOFI (Son of 
ISAAC; Moorwood et al. 19981 at the ESO’s New Tech¬ 


nology Telescope (NTT) during the nights of April 13 and 
19, 2014, under the programme number 093.C-0050(A). We 
used two low resolution grisms, covering the wavelength 
range 0.95 - 1.64 /jm (Blue Grism) and 1.5 - 2.5 /im (Red 
Grism) with 1" slit, resulting in spectral resolution R ~ 600. 

Data reduction was performed by combining the SOFI 
pipeline recipes, our IDL routines, and the IRAF task apall. 
The data reduction steps include cross-talk, flat field, and 
distortion corrections. Pairs of frames at different nodding 
positions were subtracted one from another, shifted, and 
combined into a final frame. We extracted the spectra us¬ 
ing the IRAF task apall and wavelength calibrated using 
the solution obtained by the pipeline. The spectra of tel¬ 
luric standard stars were reduced in the same way as the 
science frames. Telluric standards (spectral types B2V and 
B5V) show several prominent hydrogen lines in absorption, 
which we removed from the spectra by interpolation, prior 
to division with the black body spectrum at an appropriate 
effective temperature. This yields the response function of 
the spectrograph. Finally, we divided the science spectra by 
the corresponding response function. 

We also reduced the K-band acquisition images and 
measured magnitudes of 11.85 and 12.00 for April 13 and 
19th, respectively. These values were derived in compari¬ 
son with several 2MASS stars in the same field of view. We 
scaled the spectra so that the ratio of the integrated K-band 
fluxes corresponds to the flux ratio as determined from the 
photometry. We also scaled each of the two blue spectra 
so that they match the corresponding red spectrum in the 
overlapping wavelength range in the H-band. That way, al¬ 
though the two spectra are not calibrated in an absolute 
flux scale, they can be compared with each other across the 
entire wavelength regime. The full spectra are shown in Fig. 
[2] the H- and K-band are plotted separately in Fig. [3] 


2.3 Mid-infrared photometry 

ISO-Oph-50 has been observed multiple times with the 
IRAC instrument on-board Spitzer, mostly in the frame¬ 
work of the project YSOVAR (PI: J. Stauffer, see |Morales-| 
Calderon et al. (20111; Rebull et al. (20141). The pOph 
survey as part of YSOVAR has recently been published by 


Gunther et al. (2014). I 11 their paper, ISO-Oph-50 is found 


to be strongly variable with a reduced \ of 821.89, second 
highest in their sample, confirming the exceptional nature of 
this source. We compiled all available images for ISO-Oph- 
50 taken in the Spitzer channels 1 and 2 (3.6 and 4.5 pm) 
and measured the flux of the source. In total, we obtained 
104 epochs in IRACl and 20 epochs in IRAC2. Fluxes were 
converted to magnitudes using appropriate aperture correc¬ 
tions and zeropoints of 281 and 180 Jy for channels 1 and 2, 
respectively. The calibration was verified by comparing with 
the publicly available fluxes from the C2D Legacy project 
( [Evans et al. 2009). Typical uncertainties in the Spitzer 
fluxes are ±0.05 mag. We note that the YSOVAR lightcurve 


(see Fig. 13 in Gunther et al. (2014)) shows exactly the same 


features as ours. 

In addition, we collected 12 epochs of WISE data from 
the WISE All-Sky Single Epoch database. This includes pho¬ 
tometry in the filter W1 and W2 (3.4pm and 4.6 pm), which 
are equivalent to the Spitzer channels 1 and 2, but also 
W3 and W4 magnitudes at 12 and 22 pm, a wavelength do¬ 
main not covered anymore by Spitzer. ISO-Oph-50 is only 
detected at s/n > 10 in 8 out of 12 epochs in W3 and 5 
out of 12 in W4, in the following we only use these robust 
detections. Typical uncertainties of the WISE photometry 
are < 0.05 for W1 and W2, <0.1 for W3 and W4. Com¬ 
pared with these uncertainties, the object is significantly 
variable in all four WISE channels, with amplitudes from 
maximum to mininum of 0.27, 0.17, 0.47, 0.7mag in Wl- 
W4. The WISE colour W1 — W2 increases by 0.2 mag as the 
object gets fainter in Wl. 

The full mid-infrared lightcurve in the two short- 
wavelength bands is shown in Fig. [4] which combines Spitzer 
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Figure 3. NTT/SOFI spectra for ISO-Oph-50 taken on April 13th and 19th 2014; close-up view of the H- and K-band. Several 
photospheric absorption features are marked. 
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Figure 4. Mid-infrared photometry from Spitzer and WISE for ISO-Oph-50. The left panel shows the full lightcurve, the right panel a 
detailed view on the data from the year 2010 in the 3.4-3.6 pm band. 


and WISE data. We also list the photometry in the appendix 
in tables [A2l to IA4I 


3 THE COMPLEX VARIABILITY OF 
ISO-OPH-50 

ISO-Oph-50 has now been observed for more than 20 years. 
The observational record prior to this paper has been sum¬ 
marised by Alves de Oliveira et al. (20121, see in partic¬ 


ular their Fig. 5, but a convincing interpretation of the 
data is still missing. Our continuous monitoring in the near- 
infrared, our spectroscopy, plus the new mid-infrared data 
now adds significantly to the empirical constraints. In the 
following we will present the observational evidence in its 
entirety. 


3.1 Near-infrared lightcurve 

In Fig. [5] we plot the long-term lightcurve for ISO-Oph-50 in 
the near-infrared bands and colours, combining all literature 
data and our own new photometry. This is an updated ver¬ 


sion of Fig. 5 in Alves de Oliveira et al. (2012), containing 


datapoints from 2MASS (Skrutskie et al. 20061, UKIDSS 


1 Lawrence et al. 

2007) and a variety of other sources ( 

Bar- 

Sony et al. 1997 

Alves de Oliveira & Casaii 

2008 

Geers 

et al. 

2011 


Barsony et al. 

2012 |Scholz et al. 

2013) 

: 

The 


photometric records shows that the object is strongly vari¬ 
able. Its magnitude changes roughly from 12 to 16 in the 
K-band, from 13 to 16 in the H-band, and from 14 to 17 in 
the J-band. 


2 In addition, the object was not detected by DENIS in June 
1999, which roughly implies J > 16 and K > 14. We note that 
the object is not covered by the second release of the VISTA 
Hemisphere Survey; later releases of the survey are expected to 
add another epoch. 
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Figure 5. Long-term lightcurve for ISO-Oph-50, including data 
from various sources in the literature and our SMARTS photom¬ 
etry. See Sect. [3T] for more details. 
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Figure 6. Colour-magnitude plot for the new SMARTS near- 
infrared photometry of ISO-Oph-50 covering the years 2012-2014. 
An extinction vector is overplotted, using the reddening law by 
|Mathis| ( | 1990|) . 


There is some evidence for a bimodal behaviour, in the 
sense that the object spends prolonged periods of time near 
its brightest and faintest state. Most of our own SMARTS 
datapoints from 2013-14 are near the brightest state around 
I\ ~ 12. The object was similarly bright in 2006 in the 
lightcurve presented by [A lves de Oliveira fe Casalij ( |2008| ), 
with K ~ 12 — 13 and H ~ 14. An individual datapoint 
from 1993 is also at K ~ 12. On the other hand, the object 
was found to be significantly fainter in 2005 with K ~ 15 
and H ~ 16. Datapoints from 1999 and 2007 are on a sim¬ 
ilar level. The near-infrared lightcurves show that the tran¬ 
sition from faint to bright and vice versa has to occur on 
a timescale significantly shorter than one year. The large 
drop at the end of 2013 may indicate that the transition 
is as short as one month. In addition, Fig. [5] may indicate 
a gradual brightening between 2009 and 2014. Considering 
the years with photometry, the object seems to spend about 
half of the time in bright state and the other half in faint 
state. 

The large-scale variability evident in the long-term 
lightcurve is associated with an intriguing colour trend, first 
noted by |Alves de Oliveira &; Casali| ( |2C)0 8). The object is 
usually significantly redder when bright, or conversely, bluer 
when faint. For example, between 2005 and 2006 ISO-Oph- 
50 became brighter by several magnitudes and redder by 
0.8 mag in H — K. Also, in 2013-14 the object was bright 
( K ~ 12) and red (J — K = 3 — 4), while in 2005, 2007 
and 2008 it was faint (K ~ 15) and significantly bluer 
(J — K ~ 2). This behaviour, however, does not hold for 
all datapoints (see for example the 2MASS datapoints in 
1999 when the object was faint and red). The colour trend 
is further discussed in Sect. [5] 

Near its maximum flux, ISO-Oph-50 shows some addi¬ 
tional variability on timescales of weeks with amplitudes less 
than 1 mag in K-band, as seen in the 2006 lightcurve from 
Alves de Oliveira fe CasaII| (2008| and in our SMARTS pho¬ 


tometry from 2012-14 (see Fig. [Tj) . Our colour-magnitude 
plot (Fig. § demonstrates that during that period the ob¬ 
ject mostly becomes redder when fainter, i.e. a colour trend 


opposite the one seen in the long-term evolution. The colour 
variability is following roughly the extinction vector, which 
is overplotted in Fig. [6] This implies that variable extinc¬ 
tion is likely the main driver for the bright state fluctuations 
observed for ISO-Oph-50. In the 2013-14 data, the J — K 
colour varies by 1.4 mag, which corresponds to optical Ay 
variations by about 7 mag. Again there are two interesting 
exceptions in 2012-14 which do not fit the reddening vector. 
Compared with the reddening trend in the remaining (J- 
K,J) dataset, these two measurements are much too blue for 
the K-band brightness. This could be the same ’blue when 
faint’ effect seen in the long-term lightcurve (see above). 


3.2 Mid-infrared lightcurve 

The mid-infrared photometry from Spitzer and WISE (Fig. 
[4| shows evidence for variations by up to 3 mag at 3.6)im 
and 2 mag at 4.5 /rm, and also for variability at 12 and 22 )im 
(albeit with sparse sampling). This was already pointed out 
in jScholz et a l.||2013). Overall, the object gradually became 
brighter from 2010 to 2013, which is similar to the trend 
seen in the near-infrared lightcurve over this timespan. 

In the year 2010 the object was frequently visited by 
WISE and Spitzer at 3.6 /im, resulting in the only high- 
cadence lightcurve for this object thus far (see Fig. [4] right 
panel). The Spitzer dataset shows up to one magnitude 
of variations with a typical timescale of 1-2 weeks. These 
changes can be called ’quasi-periodic’, in the sense that the 
brightness goes up and down with some regularity, but both 
the amplitude and the period changes. This morphology 
could correspond to a series of eclipses of varying depths, 
with a decrease in flux by 0.2-0.8 mag followed by an in¬ 
crease, without flat bottom. The lack of colour information 
makes it difficult to constrain the nature of these variations 


any further. We note that Gunther et al. (2014) find a ’char¬ 


acteristic timescale’ of 5.6 d from an autocorrelation analy¬ 
sis of the 3.6 fj,m lightcurve of ISO-Oph-50, confirming that 
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the dominant timescale of this quasi-periodic behaviour is 
around one week. 


3.3 Near-infrared spectroscopy 

Our near-infrared spectra (Fig. [ 2 ]) show the typical spectral 
signature of a young mid M dwarf, with the broad peaks in 
H- and K-band caused by water absorption as well as the 


CO absorption bandheads at > 2.3 fj,m (e.g. Scholz et al. 


. The K-band shows lines at 2.21, and 2.26 /rm (see 


Fig. I3|) , which are most likely caused by Nal and Cal ab¬ 


sorption (Cushing et al. 20051. Both lines are typical for 


M3-M5 dwarfs. A few more possible photospheric lines are 
seen in the H-band (FeH at 1.62, All at 1.67, and Mgl at 
1.71 fim); all three are again characteristic for the M1-M5 
regime. 

Based on the K-band photometry in the acquisition im¬ 
ages, the object was near its brightest state when the spectra 
were taken. As pointed out in Sect. |3.1| the near-infrared 
emission in bright state is subject to significant redden¬ 
ing. To deredden the spectrum, we adopt an extinction of 
A v ~ 10, which was estimated from the K-band magnitude 
of ~ 12 in the acquisition image and the colour trend seen 
in Fig. [ 6 ] assuming a photospheric J — K of 1.0, which is 
appropriate for such an M-type object. This was done fol¬ 
lowing the typical extinction vector shown in Fig. [ 6 ] using 
the reddening law by Mathis (1990). After dereddening the 


spectra, we estimate spectral types using two different in¬ 
dices and obtain M 6 (H20 index; Allers et al.|2007 1 and M4 
(H20-K2 index; Rojas-Ayala et al. 20121. As shown recently 
by Dawson et al. (2014J7 these two indices yield consistent 
results within one subtype for young mid to late M-type ob¬ 
jects, but while the H20 index is only calibrated for spectral 
types of M5 or later, the H20-K2 index covers the full M- 
type regime. Thus, we put more credence into the result from 
H20-K2 and conclude that the spectral type is M4 with an 
uncertainty of one subtype. For comparison, the only pre¬ 
vious spectrum available was taken in the optical by |Alves| 
de Oliveira et al. ( 20121 . They estimate an extinction of 


Ay = 4.7 mag and a spectral type of M3.5, consistent with 
our value. 

Our spectra show some evidence for ongoing accretion, 
based on the presence of Paschen /? and Bracket 7 emission 
features ( Natta et al.|2004 |. This further confirms the youth 
of the object. Interestingly, the optical spectrum in faint 
state does not show any Ha emission or other evidence for 
accretion (Alves de Oliveira et al.j2012l. 

Our two spectra, taken 6 days apart near lightcurve 
maximum (K = 11.85 and 12.0), show significant spectral 
variability (Fig. [ 2 J. The spectrum taken on April 13th 2014 
is brighter and bluer than the one from April 19th. This 
variation can approximately be accounted for by variable 
extinction. Assuming Ay = 10 for April 13th and Ay = 
11.5 for April 19th would result in quite similar dereddened 
spectra with remaining differences in the range of 10 % (see 
Fig# As pointed out in Sect. |3.l] such extinction variations 
are consistent with the photometric lightcurve in the bright 
state. 


4 THE CENTRAL OBJECT: AN 

UNDERLUMINOUS VERY LOW MASS 
STAR 


With the available spectroscopy and photometry we can 
place constraints on the nature of the central object. The 
estimated spectral type of M4 would correspond to an effec¬ 
tive temperature of about 3400±200 K ( |Muzic et al.| [20l4| 
and a mass of around 0.4Mg, assuming an age of lMyr 
(BCAH isochrone; Baraffe et al. 1998). Given the uncer¬ 
tainty in spectral typing and conversion to temperature, the 
plausible range for the mass is 0.2-0 .6 Mg. Thus, the central 
object is a low mass star. Given the evidence for variability, 
accretion, strong extinction and excess emission from the 
disk, plus the uncertainty in the age and the evolutionary 
tracks, a more precise estimate of the mass is difficult to 
obtain. 

According to the 1 Myr BCAH isochrones, the K-band 
brightness for a 0.2-0.6 Mg star at the distance of pOpli 
(140 pc) and at an age of 1 Myr would be 8 to 9.5 mag. This 
is 5-7 mag brighter than the minimum light of ISO-Oph-50; 
as already pointed out by Alves de Oliveira et al. (2012) 
the object is significantly underluminous. Thus, during the 
faint state most of the light from the central source has 
to be blocked by optically thick material along the line of 
sight. A plausible explanation is that the plane of the disk 
is approximately in the line of sight, i.e. we see the object 
through an edge-on disk. The colours near minimum are only 
slightly redder than expected for a mid M dwarf (according 
to BCAH, H - K ~ 0.2 - 0.3, J - K ~ 1.0), in line with 
the presence of an edge-on disk, as in this case the radiation 
near minimum is expected to be dominated by scattered 
light, which makes the object appear bluer (for reference, 
see the models by (Kenyon et al.| 1993)). This is not unlike 
the color-magnitude trends seen in other young sources seen 
through an edge-on disk (Scholz et al.|[2010 Jayawardhana 
et al |2002 1 . 


Even near maximum the object appears underluminous 
by 2-4 mag in K and J, respectively. Extrapolating from the 
colour trend seen in Fig. [ 6 ] to a photospheric J — K = 1.0 
results in a dereddened K-band brightness of I\ ~ 10.5, 
much closer to the expected value but still underluminous. 
This could mean that the source is only partially visible in 
bright state or that is in fact at the low-mass end of the 
plausible range of masses. 


5 THE ORIGIN OF THE VARIABILITY 

5.1 Summary of observational evidence 

According to the analysis in the previous section, ISO-Oph- 
50 shows the following characteristics: 

• The central source is a very low mass star affected by 
strong reddening and excess infrared emission. 

• The object is underluminous at all times, but particu¬ 
larly so when near lightcurve minimum, by up to 5-7 mag in 
the K-band. This is indicative of an object seen through an 
edge-on disk. 

• The near-infrared lightcurve shows large-scale varia¬ 
tions over timescales of months and years. During these vari¬ 
ations the object mostly becomes bluer when fainter. 
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NTT/SOFI spectroscopy 
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Figure 7. Dereddened spectra for ISO-Oph-50 taken on April 
13th and 19th 2014. 


• Furthermore, the object exhibits small-scale changes in 
the near-infrared fluxes on timescales of weeks to months. 
Here the object becomes redder when fainter, consistent 
with changes in exctinction Ay by several magnitudes. 

• In the mid-infrared the object show variability char¬ 
acteristics similar to the near-infrarad. High-cadence mon¬ 
itoring reveals quasi-periodic changes on timescales of 1-2 
weeks. 

• There is evidence for accretion when the object is in 
maximum light, but not in minimum. 

In the following we will investigate possible scenarios that 
may explain these properties. We emphasise that what we 
propose here is mostly a qualitative scenario for the main 
features of the source; detailed modeling is required to clarify 
details. 


5.2 A clumpy disk seen at high inclination 

The most conservative statement about the cause of the 
variability is that it arises in the inner parts of the disk. 
This is supported by three pieces of evidence: a) Near- and 
mid-infrared lightcurve show the same trend, b) Short-term 
colour variations are consistent with changes in extinction 
along the line of sight; the timescales of weeks to months in¬ 
dicate that the material responsible for these changes has 
to be close, but not on the central object, c) Over long 
timescales the objects becomes redder when brighter; the 
opposite would be expected for any source of variability at 
the photospheric level (i.e. hot or cool spots). As pointed 
out in Sect. [4] the underluminosity as well as the colour 
near minimum can be understood by a disk seen at high 
inclination (i.e. close to edge-on). 

We suggest therefore that the most likely explanation 
for the complex variability of ISO-Oph-50 is a very low mass 
star surrounded by an inhomogenuous (’clumpy’) edge-on 
disk. Near its brightness maximum, the central source is 
directly visible, but seen through substantial and variable 
extinction of Ay = 10 — 15 mag. The light from the star 
may have to travel through optically thin parts of the disk, 
either the upper layers or a gap in the midplane. In addition, 


the inner disk is visible to the observer, i.e. the source also 
appears bright in the mid-infrared. In minimum, the direct 
light from the central object, the thermal emission from the 
inner disk, and the emission from the accretion flow as well 
are blocked by optically thick material along the line of sight. 
Instead, the object is seen in scattered light and thus appears 
bluer than in bright state. The change from bright to faint 
could be caused by a clump of material moving into the line 
of sight due to the disk rotation or due to precession of the 
disk. This scenario can qualitatively explain all the features 
listed in Sect. 15.11 

If the inhomogeneities in the disk of ISO-Oph-50 are 
persistent features in Keplerian rotation, the timescale of 
the variations gives us some idea about the location of the 
features relative to the central source. The quasi-periodic 
changes seen at 3.6 pm in 2010 occur with a typical cycle of 
1-2 weeks. The small-scale variations seen in 2012-14, mostly 
consistent with variable extinction, also occur on timescales 
of a few weeks. This would correspond to a feature at a 
distance of 0.05 to 0.15 AU in Keplerian rotation. For a low- 
mass star, this distance is well beyond the radius of magne- 
tospheric truncation and just beyond the dust sublimation 
radius, i.e. the features would be placed close to the inner 
edge of the dusty disk. 

However, if the variability is caused by only one feature, 
we would expect to see a repeatable pattern, stable ampli¬ 
tudes or perhaps a flat portion of the lightcurve, and not 
the irregular variability discussed in Sect. [3] More likely, a 
multitude of evolving features in the inner disk of different 
size and density is responsible for the complex variability. 
Some of them may be optically thick and thus obscure the 
central source when passing through the line of sight. Others 
are optically thin and only affect the light through variable 
extinction. With multiple features, the distance from the 
central object could be significantly larger than the value 
derived above for the single-occulter scenario. The fact that 
ISO-Oph-50 is variable all the time tells us that the in¬ 
homogeneities are present in a major portion of the inner 
disk. Furthermore, the lack of pattern indicates that some 
of the clumps may be transient features that are formed on 
timescales of weeks to years and decay on similar timescales. 

The large amplitude of the variations, the lack of a clear 
long-term period, the seemingly erratic behaviour, the vari¬ 
ous irregularities, and the colour trend observed in the long¬ 
term variations distinguish this object from other young low- 
mass stars with evidence for obscuration by disk material, 
for example KH15D ( Hamilton et al.||2001|), AA Tau (Bou 


vier et ak|2003 20131, or WL4 ( Plavchan et al.||2008 |. Fla¬ 


herty et al. (20131 find that a majority of stars with disks 


shows mid-infrared variations, albeit with much smaller am¬ 
plitudes than ISO-Oph-50, concluding that rapid changes 
in the inner disk are quite common. Based on near-infrared 


lightcurves, Carpenter et al. (2001) and, more recently, Wolk 


et al. (20131 also find objects with variability indicating 


changes in the inner disk. Some of the objects discussed 
by Wolk et al. ( |2013| ) show puzzling changes in the colour 
trend, similar to ISO-Oph-50 (see for example RWA 2 and 
RWA 7 in their paper). 

The behaviour of ISO-Oph-50 is in some ways reminis¬ 


cent of UX Ori-type Herbig Ae/Be stars (e.g. |Grinin et al. 
|1991[|Natta fe Whitney|2000| |Dullemond et al.|2003| ). These 

objects are usually described as stars with nearly edge-on 
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disks in which clumps of dust and gas are temporarily in 
the line of sight and cause irregular, long eclipses. Since the 
obscuration events in UX Ori stars occur on timescales of 
weeks to months, the occulting clumps have to be located 
in the puffed-up inner rim of the disk, rather than the outer 


flaring disk (Dullemond et al. 2003 see their Fig. 1). ISO- 
Oph-50 could be a version of the same phenomenon, but 
with a central object that is by an order of magnitude less 
massive. For such low-luminosity objects, the scale of the 
disk is different. The inner edge of the dusty disk is closer 
to the central object (<< 0.1 AU). Flaring dominates the 
disk geometry between 0.1 and 1AU (see the modeling in 
Scholz et al. J200T, [2009). Thus, in contrast to its more mas¬ 


sive counterparts, the occulters in the disk of ISO-Oph-50 
are not necessarily located at the inner edge of the disk, they 
could be spread over a range of distances. 

Thus, compared with the literature, ISO-Oph-50 is an 
extreme case of YSO variability, but following trends which 
have been observed on a smaller scale in other young sources. 
The physical origin of the inhomogeneities is difficult to con¬ 
strain from the existing data. The clumps could be density 
enhancements (e.g, spiral structure or vortices) or geometri¬ 
cal anomalies (e.g., warps or walls), as found for other young 
stellar objects. It is possible that the extreme behaviour of 
ISO-Oph-50 is solely due to a high inclination of the disk. 
Other explanations include a recent perturbation of the disk, 
for example by an accretion burst or by an instability. 

5.3 Alternative scenarios 

Here we briefly consider two possible alternative scenarios 
proposed in the literature for the behaviour of ISO-Oph- 
50. Based on the flux increase between the 2005 and 2006 
lightcurve, |Alves de Oliveira & Casalij (2008| conclude that 
ISO-Oph-50 could be an object undergoing an accretion 
burst. Since it was observed to be bright before, this would 
be a recurrent phenomenon, perhaps comparable to the pro¬ 
totype EX Lupi. Including the more recent data, there are 
several arguments against this scenario. First, as already 
pointed out in |Scholz et al.| ( |2013| ), EX Lupi is observed in 
outburst only about 10% of the time (Herbig 20071, whereas 
ISO-Oph-50 could be near maximum about half of the time 
(see Sect. |3.1[). Second, the source is bluer when fainter, 


which is not what is observed in EX Lupi-type objects (As- 


pin & Reipurtli 2009). Third, ISO-Oph-50 is clearly too faint 
for its spectral type, even in its bright state (Sect.j4|, which 
essentially excludes that the faint state represents the qui¬ 
escent period between two accretion bursts. And fourth, the 
object does not bear the typical spectroscopic signatures of 
an FU Ori burst (high veiling, no photospheric absorption 
lines, excess CO absorption and emission, see Conne fley fe| 
Greene (2010)). For all these reasons, accretion bursts alone 
do not provide a satisfactory explanation for the observed 
variability, but could contribute to some of the variations. 

Alves de Ol iveira et al.| (j2012| speculate that ISO-Oph- 
50 consists of ’two objects with similar spectral type, one 
having strong IR excess and undergoing variable circum- 
stellar extinction by the disk (i.e., close to edge-on), the 
other component either disk free or at least not obscured 
by its disk.’ At maximum brightness, both stars would be 
visible. In minimum, the primary would be obscured by its 
disk, the light thus be dominated by the secondary. Since 


the secondary is disk-less, the emission would become bluer. 
As outlined above, changes in circumstellar extinction do 
explain parts of the variability, but there is no independent 
evidence for ISO-Oph-50 being a binary. If it is a binary, we 
would expect the object to be overluminous near maximum 
when both components are visible, whereas in reality it may 
still be fainter than the intrinsic brightness of one star with 
the appropriate spectral type (see Sect. [ 4 ]). The assumption 
that the system contains a disk-less companions adds com¬ 
plications to the proposed scenario, but cannot be excluded 
based on the available data. 


5.4 Further observations 

Continued and more detailed observations of ISO-Oph-50 
could provide a wealth of information about the spatial dis¬ 
tribution of dust in the inner disk. The structures discussed 
here are too small for imaging with high spatial resolution, 
and the object too faint for optical/infrared interferome¬ 
try. The dominance of scattered light in faint state could 
be verified with polarimetric observations, but that as well 
will be challenging for this faint source. If the proposed sce¬ 
nario outlined above is correct, one might also expect signif¬ 
icant changes in spectroscopic features in the mid-infrared 
(e.g., ice feature, see Pontoppidan et al. (20051). Apart 


from these specific tests, the best strategy to map the cir¬ 
cumstellar environment of ISO-Oph-50 is prolonged multi¬ 
band (or spectrophotometric) monitoring in the near/mid- 
infrared. To test the structure of the outer disk, complemen¬ 
tary submm/mm interferometry would be beneficial. Since 
the object is currently near its brightest state, it currently 
presents us with a unique chance to obtain more detailed 
insights. 
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Table A2. WISE multi-epoch photometry for ISO-Oph-50. Only values with signal-to-noise ratio > 10 listed. 
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Table Al. SMARTS photometry for ISO-Oph-50. Typical errors 
are 0.3 mag in I, 0.1 mag in J, and 0.05 mag in H and K. 
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Table A3. Spitzer/IRAC 3.6 // photometry for ISO-Oph-50. Typ¬ 
ical errors are 0.05 mag. 


MJD 

d 

3.6 (im 
mag 

MJD 

d 

3.6 t-iTYl 
mag 

53071.3221 

13.92 

55320.5253 

12.48 

53092.2178 

13.35 

55320.8490 

12.53 

55298.4622 

11.98 

55321.2717 

12.53 

55298.7187 

12.05 

55321.7550 

12.64 

55298.9558 

12.10 

55322.3612 

12.48 

55299.3698 

12.24 

55323.0907 

12.50 

55299.8427 

12.43 

55323.7057 

12.68 

55300.3675 

12.28 

55323.9537 

12.59 

55300.8875 

12.25 

55324.2156 

12.58 

55301.6086 

12.24 

55324.4720 

12.70 

55302.1721 

12.35 

55324.8390 

12.74 

55302.2817 

12.36 

55325.2917 

12.84 

55302.5990 

12.52 

55325.8245 

13.05 

55302.8846 

12.52 

55326.4895 

13.07 

55303.2533 

12.50 

55327.1309 

12.99 

55303.7152 

12.52 

55327.2521 

12.96 

55304.2878 

12.52 

55327.4740 

12.97 

55304.9309 

12.56 

55327.7916 

13.15 

55305.4586 

12.46 

55328.1405 

13.21 

55305.5766 

12.43 

55328.5600 

13.11 

55305.8023 

12.37 

55329.2289 

13.26 

55306.1441 

12.29 

55329.8696 

13.13 

55306.6025 

12.18 

55330.4020 

13.08 

55307.1870 

12.15 

55330.6184 

13.07 

55307.7318 

12.17 

55330.9363 

12.92 

55308.4711 

12.21 

55331.2188 

12.80 

55309.1426 

12.16 

55331.5915 

12.75 

55309.2837 

12.14 

55332.0153 

12.68 

55309.6008 

12.07 

55332.5999 

12.70 

55310.0027 

12.08 

55333.3515 

16.74 

55310.4749 

12.15 

55333.9147 

12.49 

55311.0558 

12.27 

55671.0198 

11.98 

55311.6647 

12.41 

55674.5251 

12.07 

55312.2659 

12.37 

55675.5156 

12.11 

55312.9024 

12.40 

55679.0441 

12.18 

55313.0397 

12.38 

55679.9347 

12.17 

55313.3508 

12.34 
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13.33 
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55319.4430 
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56450.8949 

11.04 
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10.93 
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Table A4. Spitzer/IRAC 4.5 fj, photometry for ISO-Oph-50. Typ¬ 
ical errors are 0.05 mag. 
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13.14 

55467.3647 
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11.35 
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11.31 
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56795.8969 

10.49 
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